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We show the existence of an isolated soliton excitation over the topological ground state con- 
figuration in chiral helimagnet with the Dzyaloshinskii-Moryia exchange and the strong easy-plane 
anisotropy. The magnetic field perpendicular to the helical axis stabilizes the kink crystal state 
which plays a role of "topological protectorate" for the traveling soliton with a definite handedness. 
To find new soliton solution, we use the Backhand transformation technique. It is pointed out that 
the traveling soliton carries the magnon density and a magnetic solition transport may be realized. 



I. INTRODUCTION 

A study of spatially localized excitations over non- 
trivial many body "vacuum" configuration is one of the 
most challenging problem in condensed matter physics. 
Of particular interest is a collective transport of observ- 
able quantities accompanied with a sliding motion of in- 
commensurate (IC) phase modulation of charge and mag- 
netic degrees of freedom. However, well-known types of 
sliding density waves such as the charge-density wave 
and the collinear spin-density-wave cannot easily be ob- 
served because the internal phase modulation does not 
carry directly measured quantity^ Another example is 
the charge and spin soliton transport in conjugated poly- 
mers where the double degeneracy of ground state con- 
figurations gives rise to diffusive solitons.— In this pa- 
per, we demonstrate a new possibility to create an iso- 
lated magnetic soliton over the topological ground state 
configuration in chiral helimagnet with the antisymmet- 
ric Dzyaloshinskii-Moryia (DM) exchange and the strong 
easy-plane anisotropy. We show that the isolated mag- 
netic soliton exists and it can carry the observable mag- 
netic density. 

Hclimagnetic structures are stabilized by either frus- 
tration among exchange couplings^ or the DM relativis- 
tic exchanged In the latter case, an absence of the ro- 
toinversion symmetry in chiral crystals causes the Lif- 
shitz invariant in the Landau free energy. Consequently, 
a long-period incommensurate helimagnetic structure is 
stabilized with the definite (left-handed or right-handed) 
chirality fixed by the direction of the DM vector. Over 
the past two decades, the ground state and linear exci- 
tations of the chiral helimagnet under an applied mag- 
netic field have been a subject of extensive studies from 
both experimental and theoretical viewpoints. 5 Concern- 
ing nonlinear excitations in this class of spiral systems, 
the solitonic excitations were studied in the case of easy- 
axis anisotropy.— It has been revealed that the first sim- 
plest soliton solution in the easy-axis helimagnet presents 
a helical domain wall, a nucleation of a helical phase, 
whereas the second solution, so-called a wave of rotation, 
describes a localized change of the phase with a finite ve- 



(a) Nucleation of the spiral phase 




(b) Wave of rotation 




FIG. 1: Nonlinear excitations in a helimagnet with an easy 
axis anisotropy: the nucleation of the spiral phase (a) and the 
localized wave of rotation (b). 



locity, which is accompanied by a coming of moments out 
of the basal plane (Fig. 1). It turns out that the energy 
of the wave of rotation is less than that of spin wave with 
the same linear momentum. In a presence of magnetic 
field applied along the easy-axis, when a simple helimag- 
netic structure transforms into the conical one, there are 
again two solitons with energies less then the energy of 
spin wave with the same momentum. 

Recent progress on synthesis of new class of helimag- 
netic structures revives the interest to the case of chiral 
helimagnet with an easy-plane anisotropy. 7 Findings of 
solitons in the previous (easy-axis) case was based on 
a deep analogy between a dynamics of these nonlinear 
excitations spreading over the easy axis in the chiral he- 
limagnet and a dynamics of nonlinear excitations in the 
easy- axis ferromagnet. The correspondence reached by 
a gauge transformation enables to use a well established 
classification of nonlinear excitations in the last system^ 
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However, the theoretical tool turns out to be inappro- 
priate to be applied to the helimagnet with the Lifshitz 
invariant an the easy-plane anisotropy because the axial 
symmetry is lost in this case. 

In this paper, we show that this problem is resolved 
with the Backlund transformation (BT) technique. 10 The 
method we use has been effectively applied to studies 
of topological vortex-type singular solutions of the ellip- 
tic sine-Gordon (SG) equation^ In particular, one- and 
two-dimensional vortex lattices on both a homogeneous 
and periodic backgrounds have been constructed using 
the Backlund transformation^ One of the main findings 
of our investigation is an appearance of a non-trivial trav- 
eling soliton which carries a localized magnon density. 
This result may be useful in spintronics technology . 13 ! 14 
In Sec. II, we describe the model Hamiltonian and ba- 
sic equations. In Sec. Ill, we present the method to 
derive the novel soliton solution by using the Backlund 
transformations. In Sec. IV, we discuss the energy and 
momentum associated with a creation of the soliton. In 
Sec. V, we demonstrate that the traveling soliton carries 
the magnon density. Finally, we summarize the results 
in Sec. VI. 
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II. BASIC EQUATIONS 



We describe the chiral helimagnet by the continuum 
Hamiltonian density, 



n=-{ — \ +D[M x ^r JL -M, 



2 V dz J 

+ f3 2 M 2 + h x M x , 



dz 



dz 



(1) 



where the symmetric exchange coupling strength is given 
by a > 0, the mono-axial DM coupling strength is given 
by D, the easy-plane anisotropy strength is given by 
/3 2 > 0, and the external magnetic field applied per- 
pendicular to the chiral z axis is h = {h x ,0, 0). When 
the model Hamiltonian ([1]) is written, one implicitly 
assumes that the magnetic atoms form a cubic lattice 
and the uniform ferromagnetic coupling exists between 
the adjacent chains to stabilize the long-range order. 
Then, the Hamiltonian is interpreted as an effective one- 
dimensional model based on the interchain mean field 
picture. 



eff 



5H/5Mi 

(i = x,y, z) is the effective field acting on the magnetic 
moment M, the constant a is a feature of the symmetrical 
exchange coupling. In the explicit form these equations 



Hereinafter, we set a = 1, and rewrite h x = b, and D = 
a/2. The magnetic kink crystal phase is described by 
the stationary soliton solution minimizing TL with keep- 
ing 9 = tt/2. The solution obeys the SG equation, y zz = 
— bsimp, and is given by sin(<^/2) = sn(Vbz, q), where 
sn is the Jacobi elliptic function with the elliptic mod- 
ulus q (0 < q 2 < 1). The magnetic kink crystal phase 
is sometimes called the soliton lattice (SL) state. The 
background material behind this solution is summarized 
in Appendix A. As the modulus q increases, the lattice 
period of the kink crystal increases and finally diverges 
at q — > 1, where the incommensurate-to-commensurate 
(IC-C) phase transition occurs. 

Now, our goal is to find out possible nonlinear excita- 
tions with small deflections of spins around the basal xy 
plane. For this purpose we use the expansion 



9(z,t) = - + s6x(z,t) 



(5) 



with the small fluctuation 9\ <C 1 with 3 being a dummy 
variable controlling an order of the expansion. Plugging 
this into Eqs. (|3l4[ ). we have 



Starting with (fT]), one derives the Landau-Lifshitz an d 
equation dM. /dt = [M x H e ff], where H t 
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= 0, 
(7) 
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where the terms linear in s are hold. 

In real materials^ the order of magnitude of the DM 
coupling gives the long-period wave vector Qq ps D ~ 
10~ 2 and one obtains <p z ~ Qo ^ 10 -2 . Magnetic field 
strength b used in experiment are of order 10~ 5 . The 
further analytical treatment is performed in a regime of 
the strong easy-plane anisotropy, i.e. /3 2 > 10~ 3 . For 
spin configurations, where the last term in Eq.([7]) may 
be neglected (see the end of Sec. Ill), one obtains the 
relation 



1 dtp 
'2S/3 2 ~dt' 



(8) 



which establishes a conjugate relation between the dy- 
namical 6 and tp variables. The same relation was 
discussed in the context of soliton dynamics in one- 
dimensional magnets f§ Plugging Eq.© into Eq.©, we 
have the (l+l)-dimensional SG equation, 



bsimp 



1 d 2 p d 2 tp 



2(3 2 dt 2 dz 



= 
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We use Eqs.© and ([9|) to find the soliton solutions. 



III. BACKLUND TRANSFORMATION. 

The BT is known to be a powerful method to system- 
atically construct non-linear solution in a certain class 
of partial differential equations. The case of the sine- 
Gordon equation is briefly summarized in Appendix B. 
In the new coordinates y2b/3t — > t, \fb~z — > z , two solu- 
tions tp and tp of the SG equation tp t t — p zz = simp are 
related via the BT with the real valued parameter k^ 



<p~t — fz + k sin 
<Pz = ft + k sin 



tp+y> 
2 

tp + <p' 



1 . f tp-tp 

— sin — — 
k V 2 
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We here use the notation p z — d z tp — dip/dz, ip zz = 
d 2 <p = d 2 ip/dz 2 , and so on. A pair of the solutions are 
written in a form, 



tp{z) = 7r + 4tan" 1 [F[z)] , 
tp (z, t) = tt + 4 tan" 1 [V(z, t)} 



(11) 



It is easy to verify that the background kink crystal so- 
lution, sin((/?/2) = sn(z,q), is reproduced by choosing 
F(z) as 



F(z) = 



cn(z,q) 
1 + sn (z, q) 



(12) 



The dependence of V upon the time and the coordinate 
should be found through the BT. 

To reduce a complexity, it is convenient to rewrite 
Eqs.fTU]) in terms of the functions F and V. Then, we 



have the first and the second BTs respectively given by 



2kF z + F(k 2 + l) 2 (F 2 -l)(fc 2 -l) 
1 2k(l + F 2 ) 2k(l + F 2 ) 

2kF z - F(k 2 + 1) 
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(13) 



and 



_ V (F 2 -l)(k 2 + l) + F(V 2 -l)(k 2 -l) 
Vz ~ 2k{\ + F 2 ) • [ ' 

The further strategy is straightforward. The time depen- 
dence is firstly found from Eq. (fT3)) and then followed by 
solving of Eq. pl)) . 



A. BT Equation (Tl3|) 

The right-hand side of Eq. (fT3")) is quadratic in V. To 
reach a simplification we use the shift 



V(z,t) = U(z,t) 



(F 2 -l)(fc 2 -l) 
2 (Fk 2 + 2kF z + F) 

that transforms Eq. (TT3")) into 

U t + A(z)U 2 + B(z) =0, 

where 

Fk 2 + 2kF z + F 
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(15) 
(16) 
(17) 



and 



B(z) = [8fc(l + F 2 ) (Fk 2 + 2kF z + F) 
x [~16k 2 F 2 +F\k 2 -l) 2 
+ {k 2 - l) 2 + 2F 2 (1 + k 4 + 6k 2 )] 



(18) 



Another simplification is achieved through the identity 
(see Appendix C), 

A(z)B(z) = ±^-k 2 -±-2^=s. (19) 

In the sector s < 0, when the Backlund parameter k is 
constrained by \k + > 2/q, Eq. (|16p can be immedi- 
ately resolved 
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where s = —S 2 , and C\(z) is a function of the coordi- 
nate. The another sector s > contains no localized 
solitons. The time dependence that we need is recovered 
from Ed-ITU)]) 
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S 
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(21) 
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B. BT Equation (fl4|) 

The unknown function C\(z) should be determined 
from the second BT equation. The derivation is rele- 
gated to Appendix D and the result has the form 

A x , .. _ , . . „ , F(k 2 -l) 



A 



(At-CJ-iAj-Cu) 



(22) 



The substitution C±(z) 
equation into 
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By using the explicit expressions for F(z) and 

1 (2k \ 
A(z) = — f — dn (z, q) - (k 2 + 1) cn (z, q)\ , 

one obtain 
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After integration, this yields (see Appendix E) 



(23) 



(24) 



M{z) 



AqS 



log 



ASk~ (k 2 - l)sn(z,g) 



4Sfc + (fc 2 - 1) sn(z, q) 



4fc 2 (l-g 2 )(l + fc 2 ) 



{1 - k 2 )[(l + k 2 ) 2 q 2 



4fc 2 



k 2 - 1 

IT(n, am(z, g), q ). 

(25) 



The function IT is defined by 

fi(n, am(z,?),q 2 ) 

n(n, am(z, q), q 2 ) for n < 1, 
-U(N, am(z, g), g 2 ) + f(am(z, g), g 2 ) 
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-(l/2pi)log 



cn(z, g)dn(z, g) - p x sn(z, g) 
for n > 1, 



where am is the Jacobi elliptic function, and n = 
(1 — fc 2 ) 2 g 2 /[(l + k 2 ) 2 q 2 — 4fc 2 ] is the characteristic in- 
dex of the elliptic integral of the third kin d II. Further- 
more, N = q 2 /n, pi = y/ (n — 1)(1 — N), and is 
the elliptic integral of the first kind (see formulas 17.7.7 
and 17.7.8 in Ref. 16 ). Eqs. (|21l23l25p determine explicitly 
the solution conjugated to the soliton lattice through the 
Backlund transformation. To complete we focus on the 
asymptotic behavior of the found solution. 

For definiteness we choose |fc| > 1, the opposite case 
|fc| < 1 can be analogously treated. By noting that 
M.(z — > ±oo) = ±oo in this case, and vise versa for 
\k\ < 1, one obtains 



V(z -> ±oo, i) w T- 
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if to use the explicit expressions for F(z) and F z (z). 

Now, it is easy to prove that Eq. ppjl may be rewritten 
in the form 



V(z ~> ±oo, t) 



cn(z + 5, q) 
1 + sn(z +5,q) 



(27) 



which is similar to that used for F(z) given by Ea.(|12p. 
To find the shift 8 and explicit values of sn(<5, g), cn(<5, g), 
and dn(<5, g), we use the addition theorems for the elliptic 
functions, 



cn(z + 5) 
sn(z + 5) 



cn(z)cn(5) — sn(z)sn((5)dn(z)dn(<5) 

1 — g 2 sn 2 (z, g)sn 2 ((5, g) 
sn(z)cn((5)dn(<5) + sn(i5)dn(z)cn(z) 
1 -g 2 sn 2 (z)sn 2 ((5) : 



where the elliptic modulus g is omitted. Recalling the 
periodicity of the form (|27p and requiring the function 
given by Eq. (f2"T]) to coincide with the values of Eq. (f2"6")) 
in the boundary points and 2K of the period, we obtain 
the desired result, 



2k 



sn((5,g) =-- 

(1 + k z )q 

1 - k 2 

dn((5, g) 



(28) 



1 + k 2 ' 

From this result, we obtain the asymptotic behavior, 

cn(z + (5, g) 



(p. (z — > ±oo, t) ~ 7r + 4 tan 



1 + sn(z + 6, q) 



where the sign plus in Eqs. ([28|) is related with the limit 
z — > oo, whereas the minus is did with z — > — oo. At the 
end, we note an analogy of the introduced shift S with the 
shift of atoms relative to potential minima in the model of 
Frank and Van der Merwe (FVdM).— Now, we have done 
everything we need by using BT and found out that the 
BT creates an additional kink over the background kink 
crystal state and causes an expansion of the periodical 
spin structure at infinity. 



C. The case near the 
incommensurate-to-commensurate phase boundary 

We consider in details the case near the IC-C phase 
boundary (g — > 1). In this limit, the lattice period of the 
kink crystal tends to go to infinity and the background 
state consists of a solitary kink described by <p{z) = 
2sin~ 1 (tanhz). Fortunately, in this limit, Eq. (|24| can 
be integrated by using elementary functions. By using 
the relationships cn(z,l) = 1/coshz, sn(z,l) = tanhz 
and dn (z, 1) = 1/ coshz, Eq. ([M)l is integrated to give 



M{z) 



-z + M , 
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where A^o is a constant. The same simplifications yield 

(fc-1) 2 



A(z) 



4fc cosh z ' 



and F z (z,l) = - [cosh z(l + tanh z)] , 1) = 

[coshz(l + tanhz)] . The parameter S in Eq. ([20|) reads 
as S = (k 2 — l) / (4fc) in this case. 

After collecting the results together, one eventually ob- 
tain 



x < sinh z + cosh z tanh 



V(z,t) = 
'k 2 - 



k-1 
1 fl + k 



1 - k 



z + Mo - t 



and 



g> (z, t) = it + 4tan _1 [V(z, t)} , 



(29) 



(30) 



where v$ — (1 — k) / (1 + fc) can be thought of as the 
"velocity" of the soliton. 

The polar angle computed via Eq.© takes the form, 




(31) 



This solution is viewed as a "collision" of two kinks as 
shown in Fig 2. One of the kinks is "at rest" as a member 
of the background configuration, while the other travels 
with the speed vq and go through the background with- 
out changing its shape. 

To justify the inequality \d 2 9 1 /dz 2 \ < 2/3 2 6>i that have 
been assumed to obtain Eq.© the restriction 



< 4/3 2 



(32) 



must be fulfilled, which is derived by means of Eg. (|3Tj) . 
To obtain this condition, we use the asymptotic of 
\(d 2 9i/dz 2 )/9\\ at z — > ±oo. Note that the constraint 
\k + l/k\ > 2/q excludes only the points k = ±1. 
Therefore, there is a region of finite k values provided 
(3 2 > 1/2. 



IV. ENERGY AND MOMENTUM 

Next we consider the energy and momentum associated 
with creation of the found soliton. The energy density 
([1]) is written in polar coordinates as 



H=^(0 2 z +sm 2 e<pl) + ^sm 2 



~+P cos 0+b sin 9 cost/?. 

(33) 




FIG. 2: A profile of the soliton excitation at q — *• 1 in a spiral 
with an easy plane anisotropy directly before (a) and after 
(b) "collision" of the kinks. 



Using the expansion (|5I8|) and neglecting terms higher 
than second-order derivatives one obtains 



H 



1 



:<Pz 



4/3 2 



&COS!^ 1 



8f3 4 



(34) 



In the case of small fields, the term proportional to bip 2 
may be ignored and the energy density TL measured in b 
units becomes 



H 



where c = a/ and the coordinates V%bf3t — * t, 

\fb~z — > z are again used. Following the method sug- 
gested in Ref 3, we find the difference between the en- 
ergy densities calculated on the solutions coupled by the 
BT 



+ 



H{<p)-H{y) = c{(p z - ip z ) 
k sin " 
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k sin 
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— — sm 

k 
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ip-g> 



ft 
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where the transformations (|10p are employed. This 
amounts to a derivative of the function 



2k cos 



v + (p 



$( e ) = c {ip - ip) + % cos (- — ^ 
k \ 2 

(36) 

with respect to z. Thus, the difference of the energies 
integrated over the total length L of the system is 
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equal to 



Svq9i has the explicit form 



dz {H(<p)-H(<p)} = * (e) (L) - * (e) (0) . 



The momentum density V = hS (1 — cos#) ip z is 
treated in the same manner. The spin value S is re- 
alted with the magnetic moment by M = 2/iqS. Using 
again the expansions (|5l8p and the new coordinates one 
obtains 



of the same accuracy as in the case (|34l) . Here, the mo- 
mentum V is measured in the units hS\/b, and v = 
y / b/2f3 2 . The difference between the momentum den- 
sities of the solutions conjugated by the BT amounts to 

v(<p)-v{<p) = &r\ 

where the function vj/( m ) is given by 



ip—ip + v 



2 {<p-<p 
— cos — - — 
k \ 2 



2k cos 



ip + ip 



(37) 

Therefore, the additional momentum with reference to 
the background kink crystal state becomes 



J dz {V{(p)-V{tp)} = * (m) (L) - * (m) (0) . 



V. MAGNON CURRENT. 

The magnon current transferred by the ^-fluctuations 
is determined through the definition of the accumu- 
lated magnon density p s in the total magnon density 
Af = gp^S (1 — cos 6) — p + p s , where the "superfluid" 
part p s = — §pbS cos9 is conjugated with the magnon 
time-even current carried by the ^-fluctuations. Then, 
we obtain the magnon current via the continuity equa- 
tion Aft + J% = 0. Here, we compute the magnon current 
density carried by the traveling soliton in the limit of 
q — ► 1, where the analytical solution with the intrinsic 
boost transformation is available based on Eqs. ([30"|) and 

(ED- 

The additional magnon density associated with the so- 
lutions given by Eqs. (f50|) and (f3"Tj) is 



SAf ■■ 



S 

W 2 



Taking the time derivative of SAf and using the property 
(pt = — vq (p z we get 5 Aft = — J*, where the current J z = 



J z (z,t) = S< 



I b (1-fc 2 ) 1 
2/3 2 k 1 + V 2 
coshz 



cosh 



k A - 1 fl + k 



4fc V 1 - k 



z + Mo-t 



(38) 



Note a significant difference of the result with the magnon 
current due to the translational motion of the whole kink 
crystal considered by some of the authors recently^ In 
the latter case, the sliding motion of the whole kink crys- 
tal excites massive spin wave excitations of the 0-mode 
above the traveling state that is responsible for a magnon 
density transport. In the present case, the nontrivial soli- 
ton solution itself carries a localized magnon density (the 
magnon " droplet" ) due to the intrinsic boost symmetry. 

In Fig. 3(a), we depict the background topological 
charge Q = d z <p associated with the standing kink around 
z = in the q — > 1 limit, i.e., ip(z) = 2sin _1 (tanh z). In 
Figs. 3(b-l)-(b-6), we show the magnon density distri- 
bution J z (z, t) carried by the traveling kink at the time 
t = —20, —4, —2, 0, 2, 4,, respectively, where the soliton 
travels from left to right. The traveling soliton collides 
with the standing kink at t = 0. It is clearly seen that 
the magnon density is largely amplified when the soliton 
"surfs" over the standing kink. 



VI. CONCLUDING REMARKS 

In summary, by using the Backlund transformation 
technique we investigated soliton excitations in the 
chiral helimagnetic structure with the antisymmetric 
Dzyaloshinskii-Moryia exchange and with the strong 
easy-plane anisotropy, which is experienced by the ex- 
ternal magnetic field applied perpendicular to the mod- 
ulation axis. The soliton we found was obtained as an 
output of the BT from the kink crystal solution as an 
input. An essential point is that the traveling soliton 
cannot exist without the kink crystal (soliton lattice) as 
a topological background configuration. We may say that 
the nontrivial topological object is excited over the topo- 
logical vacuum. The standing kink crystal enables the 
new soliton to emerge and transport the magnon den- 
sity. As compared with the motion of the whole kink 
crystal with a heavy mass^ our new soliton is a well 
localized object with a light mass. This object should 
be certainly, more easily triggered off and propagate over 
the crystal. 

We stress that our soliton has definite chirality, because 
of the presence of the DM term Ddip/dz in the Hamilto- 
nian density Eq. |T]) . The presence of this term lifts the 
degeneracy between the left-handed soliton and the right- 
handed antisoliton solutions. For example, in Eq. tfTTj) . 
the right-handed antisoliton solutions may be given by 
changing the sign of the phase gradient, i.e., (p (z) = 
7T - 4tan _1 [F(z)], and (p\z,t) = tt - 4tan~ : [V{z,t)}. 
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FIG. 3: (a) The background topological charge Q — d z *p 
associated with the standing kink. (b-l)-(b-6) The magnon 
density distribution J z (z,t) carried by the traveling kink at 
the time t = -20, —4, —2, 0, 2, 4, respectively. The soli- 
ton travels from left to right. The magnon density is largely 
amplified when the soliton "surfs" over the standing kink. 
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APPENDIX A: FORMATION OF THE KINK 
CRYSTAL STATE 

We here discuss the kink crystal (soliton lattice) for- 
mation from a general view point. Let us consider a 
magnetic system described by the two-component or- 
der parameter (OP) (77, £) with the Ginzburg-Landau 
functional £ 



$ = i / dz [r(r)0 + u{rtf + w(ri n + £ n ) 



+ia 



<2£ _£*7 
dz dz 



drj dt; 

^ dz dz 



(Al) 



where the condition u > 0, 7 > ensures a stability of ex- 
tremal points of the functional. The signs of w and a are 
arbitrary. We assume the spin arrangement is uniform 
in the x and y directions and hence the volume integral 
is implicitly reduced to one-dimensional integration over 
z-axis, where L is a crystal size in this direction. In the 
approximation of constant OP modulus, p = const, when 
i] = pe %v and £ = pe~ t(p , the functional (|Aip depends 
only on the phase ip, 



Although these solutions satisfy the same SG equation 
as Eq. their static energies are higher than the left- 
handed soliton solution given by Eq. lfTTj) . 

Finally, it would be of interest to discuss a difference 
between the conventional single Bloch wall and our soli- 
ton. The Bloch wall is formed within a non-topological 
background, just a ferromagnet, and it dephases easily. 
On the other hand, our soliton should be more robust 
against a dcphasing, because it emerges within the topo- 
logical background (soliton lattice). Chirality and topol- 
ogy support a stability of the moving soliton. This con- 
sequence is quite obvious in the context of the soliton 
theory, but may serve a quite new strategy in the field of 
spintronics. For example, in the left-handed chiral crys- 
tal, only the left-handed kink crystal would be formed 
and our soliton inherits with the corresponding chirality. 
Essentially, the crystallographic chirality plays a role of 
protectorate for the background chiral spin texture and 
causes the traveling soliton over the background. This 
new traveling soliton can be regarded as a promising can- 
didate to transport magnetic information by using chiral 
helimagnet. 



$ = rp 2 + up 4 



L ldZ 



IP 



^) + 2(jp 2 ^- + 2wp n cos(ntp) 
dz 1 dz 



(A2) 



and includes p as a parameter. Minimization of $ with 
respect to ip results in the equation 



if 

dz 



-j-s (ntp) + v sin(rn^) = 0, 



(A3) 



where the effective anisotropy parameter is defined by 
v = n 2 {w / p( n ~ 2 ^ . The case of magnetic field corre- 
sponds to n = 1. 

Without of the nonlinear anisotropy term, Eq. (|A3|) is 
resolved by ip — Qz which describes an one-harmonic 
IC structure, for example, a simple helimagnet, with the 
wave vector Q = — er/7. At finite v, the exact periodic 
solution is given by 



sin 



n ' 



-z,q 



(A4) 
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The elliptic modulus q must be determined by minimizing 
the corresponding energy, 



Two successive BTs commute, i.e. B ai B a . 2 = B a2 B ai , if 
the Bianchi's identity 



2p 2 7 - 



up 4 - 2p z \a 

2 



nqK 
4 E 



(A5) 



where K and E denote the elliptic integrals of the first 
and second kind, respectively. This procedure yields q as 
a function of the anisotropy parameter v, 



E/q = \/v c /v, 



(A6) 



with the critical anisotropy parameter being defined by 
v c = n 2 7r 2 cr 2 /167 2 . A change of q from to 1 corresponds 
to a change of v from to v c . Varying the parameter v 
causes a drastic change in the behavior of the amplitude 
(|A4[) . The region of an almost constant phase within 
the period I comes up at v — > v Cl while the phase rapidly 
changes at the ends of the period, where the overall phase 
change is 2-njn. The region of the constant phase in- 
creases as v — > v c . For < v < v c , the kink crystal phase 
is stabilized, where a periodic array of C-phase regions 
separated by the kinks (solitons). The spatial period is 
given by I = AqK/^/v, and it diverges logarithmically at 
v — > v c , i.e. q—*l, 



I = (4g/VU) ln^/Vl-g 2 



(A7) 



APPENDIX B: BACKLUND TRANSFORMATION 

An existence of exact multi-soliton solutions is a pe- 
culiar property of the SG equation, and the BT is a sys- 
tematic way to obtain them. Indeed, let both tpo and ipi 
are solutions of the SG equation 

d+d-ip = sin ip 

written via light-cone coordinates x + = (x + t) /2, and 
x~ = (x — t) /2. Then, the Backhand transformation 
ipi = B a [ifo] is given by 



e sin I | , (Bl) 



where a = e x is called a scale parameter. The relation is 
consistent with the SG equation, i.e., d-d+ipo — sin</?o 
and d—d+<pi = sintpi. Any two functions ipo and tpi 
that satisfy the BT necessarily solve the SG equation. 
Eq. ljBTj) is nothing but Eq. (fTU|) . 

It turns out that analytical expression for multi-soliton 
solutions may be outlined by an entirely algebraic proce- 
dure because the BT embodies a nonlinear superposition 
principle known as Bianchi's permutation theorem. Sup- 
pose that ipo is a seed SG solution, and cp\ t 2 are the BTs 
of ip a , i.e. 

<pi = B ai [ip ] , ip 2 = B a2 [ip ] ■ 



<^3 = + 4 tan 



a 2 + ai / ip 2 - <pi 

tan 

a 2 - di 



is fulfilled. It means that the non-linear superposition 
rule holds ip 3 = B a2 [tpi], ^3 = B al [ip 2 }. This alge- 
braic relation indicates that a series of soliton solutions 
is given by if = 4tan _1 [f/g], which supports the forms 
of KqJ[ID. 



APPENDIX C: COMPUTATION OF THE 
PRODUCT A(z)B(z) 

Our objective is to find the product A(z)B(z). By 
using Eqs. (|17ll8[) we obtain 

A(z)B(z) 
1 



16/c 2 (l + F 2 ) 2 

[l6fc 2 F 2 - F\k 2 - l) 2 - (k 2 - If - 2F 2 (1 + k 4 + 6k 2 )] 

1 



16k 2 (l + F 2 ) 2 
16k 2 F 2 - 4F 2 (k 2 + l) 2 - (F 2 - l) 2 (fc 2 - 1) : 



The function F{z) has the derivative 



(CI) 



F z {z) = ---. MZ,<1) 



q 1 + sn(z, q) ' 



that yields for the numerator in Eq. (|Cl 



l6k 2 F 2 - 4F 2 (k 2 + l) 2 - (F 2 - l) 2 (fc 2 - l) 2 



= 16 



fc 2 dn 2 (z,g) 
q 2 (l + sn(*,g)) 



^-4(fc 2 + l) 2 Cn2(z ' g) 



(1 + sn(z, q)) 



, 4(fc2 _ 1)2 **{z iq ) 
(l + sn(z,q)) 



(l + sn(z,q)y \ q 



4fc 2 



- fc 4 - 1 - 2fc' 



(C2) 



Thus, we reproduce the result (fl9|) . Note that the prod- 
uct A(z)B(z) embodies no coordinate dependence, hence 
it equals a constant s. 
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APPENDIX D: DERIVATION OF EQ. ([22j) 



Eq.flDT 



To derive the determining equation for the function 
C\{z) we rewrite Eq. (fl4|) through the function U (z,t) 



U z 



FF z (k 2 - 1) 



Fk 2 + 2F z k + F 
(k 2 - 1)(F 2 - 1) (F z k 2 + 2F zz k + F z ) 

2 (Fk 2 + 2F z k + F) 2 
F 2 (k 2 + 1) ( (fc 2 - l)(F 2 - 1) 



2k(l + F 2 ) 

(fc 2 + 1) 
2fc(l + F 2 ) 

F(k 2 - 1) 



U 



u 



2 (Fk 2 + 2kF z + F) 

(k 2 ~l)(F 2 -l) 
2 (Fk 2 + 2kF z + F) 



2fc(l + F 2 ) 



n (k 2 -l)(F 2 -l) 

2 (Fk 2 + 2kF z + F) 



and use Eq. (|2"0|) to obtain the relation, 



S 
A 



F 2 (k 2 + 1) A z 



(k 2 + 1) 



2fc(l + F 2 ) A 2k(l + F 2 ] 



F(k 2 -1) (k 2 -l)(F 2 -l) 



2k(l + F 2 ) (Fk 2 + 2kF z + F) 
2 1 



tanh(T) 



A J cosh 2 (T) 

(fc 2 + l)(fc 2 - 1)(F 2 - l) 2 
_ 4fc(l + F 2 ) (Fk 2 + 2kF z + F) 
(fc 2 - 1)FF Z (k 2 - 1)F 



Fk 2 + 2kF z + F 2fc(l + F 2 ) 
'S\ 2 (k 2 -l)F 
.A. 



2fc(l + F 2 ) 
(fc 2 - 1)(F 2 ^ l) ( Fz k 2 + 2kF zz + F z 
2 (Fk 2 + 2kF z + Ff 
(fc 2 - lfF(F 2 - l) 2 



8fc(F 2 



1) (Fk 2 - 
-0, 



2fcF z + FY 



(Dl) 



where T = (S/A) (At — C\). Being constant at any 
time moment, Eq. ljDip means that the coefficients before 
tanh(T) and cosh -2 (T) turn simultaneously into zero. 
With some tedious but a straightforward algebra (see 
below) one finds that only the factor before cosh -2 (T) 
produces the non-trivial result 



S 
A 



F 2 (k 2 + 1) A± (fc 2 + 1) 
2fc(l + F 2 ) + ~A ~ 2fc(l + F 2 ) 
F(fc 2 -1) (fc 2 -l)(F 2 -l) 



2fc(l + F 2 ) (Fk 2 + 2kF z + F) 

S 

" 2fc^l 2 (l + F 2 ) (Fk 2 + 2kF z + F) 

x [^(fc 2 -!) 2 ^ 2 -!)^ 

- 2fcA 2 (l + F 2 ) (F (fc 2 + 1) + 2kF z ) 

+A(k 2 + 1)(1 - F 2 ) (F (fc 2 + 1) + 2kF z )} 

S 

~ A 2 (I + F 2 ) (Fk 2 + 2kF z + F) 
x [A(l - F 2 ) (2kF + F z (k 2 + 1)) 
-A z (l + F 2 ) (F (fc 2 + 1) + 2fcF z )] . 

Now we use Eq. fTTf to find 

A z 1 1 



(D2) 



(D3) 



A 1 + F 2 F(k 2 + \) + 2kF z 
x [2fci^(l + F 2 ) + (fc 2 + 1) F z (l + F 2 
-AkFF 2 - 2F 2 F Z (fc 2 + 1)] . 

The function F(z) obeys the equation 

(1 + F 2 )F ZZ - 2FF 2 = -F(F 2 - 1), 



which is derived from the SG equation Together with 
Eq. (|D3[) this produces 

A, _ l-F 2 F z (fc 2 + 1) + 2kF 
~A ~ 1 + F 2 F(fc 2 + l) + 2fcF 2 ' 



Plugging this into the right-hand side of Eq. (|D2|) we come 
to zero. 

After lengthy manipulation it may be shown that the 
free terms in Eq. (|Dl[) . i.e. those without either tanh(T) 
or cosh -2 (T) factors, are simplified to 

- (fc 2 - 1) [g 2 fc 4 - 2fcV(8S 2 - 1) - 4fc 2 + q 2 ] 



x cn (z, q) 



[(fc 2 + 1) qcn (z,q) - 2fcdn (z, q)] ' 
16fc 3 q 5 [1 + sn( Zl q)f 



The factor g 2 fc 4 - 2k 2 q 2 (8S 2 - 1) - 4fc 2 + q 2 equals to 
zero bearing in mind Eq. (|19p . 



APPENDIX E: DERIVATION OF EQ. (p5|l 

To obtain M.(z), the derivative ((24)) is split into two 
parts, 



Indeed, consider the coefficient before tanh(T) in 



M z (z)=M (z)+M 1 (z), 
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where 
M (z) 
Mi(z) 



(1 — k 4 )q 2 cn 2 (z, q) 



(1 + fc 2 ) V - 4fc 2 - (1 - k 2 ) 2 q 2 sn 2 (z, q) ' 

2kq(l — k 2 )cn(z, q)dn(z, q) 
(1 + fc 2 ) V - 4fc 2 - (1 - k 2 ) 2 q 2 sn 2 (z, q) ' 



and both terms are separately considered. 

The integration of M.q{z) is straightforwardly per- 
formed and one obtains 



/ 



dzMo(z) 
k 2 + l 



-z + 



Ak 2 {l~q 2 )(l + k 2 ) 



k 2 -l (l-fc 2 )[(l + fc 2 )V-4fc 2 ] 



x n 

where 
U(u,a,q 2 ) 



; ,am(z,q),q : 



(i-fc 2 )V 

1 + fc 2 ) V - 4k 2 

u q 2 sn(a 1 q)cn(a, q)dn(a, g)sn 2 (w, q) 



I 



1 — (7 2 sn 2 (a, g)sn 2 (u, q) 



du 



is the elliptical integral of the third kind. 

Taking into account the relationships (1 + k 2 ) 2 q 2 — 
4fc 2 = 16q 2 k 2 S 2 and cn(z, g)dn(z, q) = (d/dz)sn(z, q) the 
second term can be written as follows 



Mi(z) 



45 



ASqk — (1 — k 2 )qsn(z, q) 



ASqk + (1 — fc 2 )gsn(z, q) 



dz 



sn(z,q) 



that yields the desired result 



/ 



dzMi(z) = ^log 



ASk- (fc 2 - l)sn(z,g) 



ASk+ (fc 2 - l)su{z,q) 
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